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Persistent currents with non-quantized angular momentum 
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We analyze the generation of persistent currents in Bose-Einstein condensates of ultracold gases 
confined in a ring. This phenomenon has been recently investigated in an experiment [Nature 506, 
200 (2014)], where hysteresis loops have been observed in the activation of quantized persistent 
currents by rotating weak links. In this work, we demonstrate the existence of 3D stationary 
currents with non-quantized angular momentum. They are generated by families of solitary waves 
that show a continuous variation in the angular momentum, and provide a bridge between different 
winding numbers. We show that the size of hysteresis loops is determined by the range of existence 
within the weak link region of solitary waves which configure the energy barrier preventing phase 
slips. The barrier vanishes when the critical rotation leads winding numbers and solitonic states to a 
matching configuration. At this point, Landau and Feynman criteria for phase slips meet: the fluid 
flow reaches the local speed of sound, and stationary vortex lines (which are the building blocks of 
solitons) can be excited inside the system. 

PACS numbers: 03.75.Lm,67.85.De,67.85.Lm,05.30.Jp 


I. INTRODUCTION 

Persistent currents are one of the striking phenomena 
associated to superfluidity. In charged systems, electric 
currents flowing without dissipation along superconduc¬ 
tor loops are due to the superfluid flow of Bose-condensed 
Cooper pairs, and have found an extensive field of ap¬ 
plications since the prediction of SQUIDs [T. In neutral 
systems, persistent currents have been observed since the 
early experiments with superfluid liquid helium (see Ref. 
[21 and references therein), and have experimented a re¬ 
vival with the advent of Bose-Einstein condensation of 
ultracold gases. Multiply connected geometries have al¬ 
lowed the generation of long-life currents in Bose-Einstein 
condensates (BECs) [3HS , paving the way for atomtronic 
devices. 

The effect of a rotating external potential producing a 
depletion in the density (weak link) has been theoretically 
and experimentally addressed in rings of ultracold atomic 
gases HEHE]. The rotating weak link drags the super¬ 
fluid into a moving state once the resistance of the system 
to generate excitations has been overcome. Whether su¬ 
perfluid excitations are triggered by reaching the speed 
of sound, according to Landau criterion, or by vortex 
nucleation, following Feynman criterion, has yet to be 
elucidated [7j. Experimental works (see Ref. [9] and ref¬ 
erences therein) point out that vortices are responsible of 
the activation of persistent currents in a ring. Phase slips 
seem to be originated by the transit of a vortices towards 
the inner (outer) edge of the system [TP] , producing the 
increase (decrease) up to one unit in the winding number. 

In Ref. [9 j a hysteresis loop has been observed in the 
generation of quantized persistent currents by means of 
a rotating weak link. A 3D ring of Bose-condensed 23 Na 
atoms was prepared in a stationary state with winding 
number 0 (1), and the frequency U+ (U~) at which the 
winding number changes to 1 (0) was measured. The 


size of the hysteresis loop, i.e. the difference (U+ — fi“), 
was compared with numerical results obtained by solv¬ 
ing the 3D Gross-Pitaevskii equation (GPE), showing a 
big discrepancy. As the intensity of the link decreases, 
the size of the loop becomes considerably smaller in the 
experiment than in the numerics. The reason adduced 
for the difference points to a dissipation mechanism, but 
a complete explanation has still not been given. 

Hysteresis is a phenomenon intimately linked to sys¬ 
tems presenting multiple local energy minima separated 
by energy barriers that must be overcome when driving 
the system from one minimum to the other m- The fi¬ 
nal state depends on the path followed in the dynamical 
evolution searching for equilibrium states. For superflu¬ 
ids in a ring, it is well known that local energy minima 
correspond to states with different winding number [42]. 
However, neither the exact nature of the states lying on 
the barrier nor their stability have been shown yet, al¬ 
though it has been hypothesized that unstable states sep¬ 
arate adjacent winding numbers. 

In this work we demonstrate the existence of station¬ 
ary states with non-quantized angular momentum in 3D 
rings. These states are solitary waves that fill the gap of 
angular momentum between winding numbers. Among 
such solitonic states, there exist dynamically stable con¬ 
figurations that can support long-life currents. In addi¬ 
tion, we show how the rotating weak link provides com¬ 
mon features of solitary waves to winding number states. 
In doing so, the size of hysteresis loops can be deter¬ 
mined by the range of existence of solitary waves within 
the weak link region. Solitonic states configure the en¬ 
ergy barrier preventing phase slips and meet a family of 
winding number states when the barrier vanishes. This 
fact indicates that Landau and Feynman criteria may 
also match. The former stablishes that phase slips are 
produced when the fluid flow reaches the local speed of 
sound, at the end of a winding number family of station- 
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ary solutions. But this is also in agreement with Feynman 
criterion, which states that phase slips appear when vor¬ 
tex lines can be excited inside the system, as it occurs at 
the end of a solitonic family of steady states. 


II. MODEL 

We will consider BECs in a ring-shaped geometry 
(torus in 3D) of radius R rotating at frequency D within 
the mean-field regime. Stationary states T(r,£) = 
fulfill the time-independent GPE (in the co¬ 
moving frame) with chemical potential g: 

(-£ v2 -^ + nr) + ^)^^, (1) 

where g = Airffia/m defines the interaction strength 
as a function of the s-wave scattering length a, and 
the external potential E(r) includes harmonic trapping 
^trap( r ) = m<jj 2 z 2 /2 + muj\(r±_ — R) 2 / 2, with r\ = x 2 + 
y 2 , and a weak link U w i( r) = Uo exp[—2 (R(9 — Oq)/wq) 2 ] 
along the azimuthal coordinate 0 = arctan (y/x). 

The angular momentum operator is defined by L z = 
— ihde , and the energy can be calculated by E[ip\ = 
f dr {|(—ifrV — mtt x r)'0| 2 /2m + V\ip\ 2 + g\ijj\ A /2}. 

Persistent currents can be supported by states with 
integer winding number g, quantized angular momentum 
per particle ( L z )/N = hq , and constant density along 6 , 
that is 


*,(r ,t)=e- i ^ q (z,r ± )e i i e . (2) 

In addition, persistent currents can also be provided by 
solitonic states, which present non-quantized angular mo¬ 
mentum and produce density dips along the ring 

^/t(r, t) m e iK$ , (3) 

where /c, a real number, plays the role of a quasimomen¬ 
tum (see below). 


III. ID RINGS 

This framework can be easily visualized in ID rings. 
Without weak link, the lowest-energy states for a given 
angular momentum, or yrast states, have been identified 
as solitary waves connecting states with different wind¬ 
ing numbers [IMS]. Assuming that a tight transverse 
harmonic trap enforces a ID geometry, and using units 
of Er = h 2 /mR 2 , the dimensionless ID version of Eq. 
0 without external potential reads 

I(-ido-fi) ip + g\ip\ 2 tlj = flip, (4) 

where jl = h/Er + D 2 /2, Cl = HQ/Er , g = 2 aRN/a 2 ± , 
and a± = h/mu:±_ is the characteristic length of the 




FIG. 1. Soliton connection between winding-number states in 
a ID ring for condensates with g = 100. Only winding num¬ 
bers 0 and 1 are considered, without (a) and with (b) external 
weak link of depth U w i = 0.56 g and width wq = R along the 
angular coordinate. For each case, angular momentum per 
particle is shown above and energy per particle below. The 
thick dotted lines of the top graph (a) indicate the values of 
the speed of sound, Co and ci, relative to the winding number 
state 0 (right) and 1 (left). The hysteresis loop is displayed 
in (b) by thick gray lines and arrows. 


transverse trap. Among the analytical solutions of Eq. 
©. a dark soliton moving with linear velocity v s along 
the ring is given by Eq. 0 with 





i/3 + ytanh 
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where y 2 + /3 2 = 1, (3 = v s /c, c = y/ g/m is the speed of 
sound, and hi^/mR = flR — v s . As can be seen in Fig. 
[lja), solitons (open circles), given by Eqs. (|3| and ([5]), 
connect states ip q (0) oc e tq0 having consecutive winding 
phases (solid and dashed lines). 

By producing localized density depletions in the ring, 
a single soliton induces a phase jump AS = 2 arccos /? in 
the condensate wave function m- These phase jumps 
need to be counterbalanced by background currents 
around the torus, with angular speed k = —q — AS/2ir, to 
ensure a single-valued order parameter whenever a closed 
loop is followed. The resulting flow carries angular mo¬ 
mentum per particle, (L z )/Nh = — — K+(Aq — 1)/2, 

which takes continuous values in the range between inte¬ 
ger winding numbers q and q + Ag, with A q = ±1. The 
velocity of moving solitons v s cannot exceed the speed of 
sound c, at which nonlinear soliton excitations transform 
into linear sound waves traveling on a winding number 
state. Therefore this speed limit bounds the range of 
existence of solitonic solutions. Provided that solitonic 
states can show dynamical stability na, their associated 
flows will appear as persistent currents. We will see that 
there exist the counterpart of these states in multidimen¬ 
sional systems. 


A. Weak link 

The presence of a rotating weak link introduces varia¬ 
tions in the density profile that give rise to phase jumps, 
and, like in solitons, background currents are generated 
in order to compensate such jumps. Fig. ljb) dis¬ 
plays our results obtained from the numerical solution 
of Eq. © for a system with g = 100, and weak link 
of depth U w i = 0.56 g and width = R. As can 
be seen, and contrary to the case without weak link, 
both winding number and solitonic states show varia¬ 
tions in the angular momentum with the rotation rate. 
This effect can be deduced from the expression for the 
current density j(r, t) = p(h\7S/m — x r), expressed 
in terms of densi ty p a nd phase S of the order param¬ 
eter T(r,£) = y/ p(r, t)e zS ^ r,t \ Defining the velocity as 
v = j/p we get v{6) = h{deS{0) — Cl)/mR for the ID case. 
Solving for the phase gradient and integrating along the 
whole ring we obtain 

27 rq = 27 rCl + j) d0v{6) , (6) 

where, due to the single-valuedness of the order parame¬ 
ter, we used S{6 = 27t) — S{6 = 0) = 2irq. Whenever the 
dimensionless angular rotation Cl takes integer values, it 
is possible to look for static states having v(6) = 0. In¬ 
stead, stirring the system at non-integer Cl values imposes 
a non-vanishing velocity, hence phase gradients. 

As a consequence of the weak link, which breaks the 
rotational symmetry of the system, a gap is opened in the 
energy spectrum of winding numbers D-p. Our numerical 
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FIG. 2. Phases (left column) and dimensionless azimuthal 
densities (right column), obtained from numerical solutions 
of the time-independent ID GPE Q, corresponding to wind¬ 
ing number q — 0 (above) and solitonic states (below) for the 
same parameters as in Fig. B b )’ rotating at angular frequen¬ 
cies Cl — 0.5, 1.0 and 1.4. Thick arrows (on the left column) 
indicate the sequence of both type of states for increasing Q, 
approaching convergence near = 1.4. 


results, displayed in Fig. 0b), quantify the qualitative 
analysis of this fact made previously in the literature. As 
can be seen, swallow tails are developed in the presence 
of a weak link, cutting the otherwise continuous curves 
traced by winding numbers in Fig. [lja). Generating 
solitons at the weak link position is energetically more 
favorable than at higher-density regions. For weak link 
extents wq greater than the healing length £ = y/h 2 /mg , 
a local speed of sound c(0) = yfgn\(6)/m can be defined 
(n i being the azimuthal density), which takes the low¬ 
est value at the weak link center. As the speed of sound 
is a measure of the range of existence of solitonic solu¬ 
tions, this range is reduced inside the low density region, 
relative to the system without weak link. 

When rotation is introduced, opposite effects can be 
observed in solitonic and winding-phase states (see Fig. 
|2j. For increasing angular frequency, solitons decrease 
their angular momentum, whereas winding numbers in¬ 
crease it. In this way, both type of solutions converge 
at some intermediate value, where the superfluid flow 
reaches the speed of sound at the weak link. Therefore, 
the bounds for the existence of solitons are also applica¬ 
ble to states with a winding phase. 

Hysteresis loops, as those observed in the experiment 
|9j, can be clearly identified in Fig. [ljb) (thick brown 
lines and arrows). Distinct from previous analyses, our 
results are able to measure the size of the energy barrier 
preventing phase slips. Solitonic states are the responsi¬ 
ble of such barrier, occupying stationary points in the en¬ 
ergy landscape of the system. The barrier vanishes when 























4 


both type of states connect, and it is no longer possible to 
continue beyond these points by increasing the angular 
velocity. As a result the system decays to the next lower- 
energy winding-phase state. When the intensity of the 
weak link grows, the speed of sound and hysteresis cy¬ 
cle width reduce, and solitons and winding-phase states 
become energetically closer. 


IV. 3D RINGS 

In multidimensional systems, the ID picture outlined 
before is reproduced: solitons are the stationary states 
making the connection between winding numbers and 
building the energy barrier separating them. In addi¬ 
tion, dimensionality allows the emergence of novel soli- 
tonic states and introduces new issues concerning their 
stability. It was demonstrated in Ref. jlT] that differ¬ 
ent transverse standing waves can be excited in straight- 
channeled superfluids. A dark soliton possesses the sim¬ 
plest structure among such standing waves, and it is a 
dynamically stable configuration for small values of the 
chemical potential, when no other transverse states can 
be excited. But more complex configurations containing 
vortex lines are also possible. They emerge at higher val¬ 
ues of the chemical potential, and thus higher densities, 
and cause the dark soliton to turn into a dynamically 
unstable configuration. The stable state is then the soli- 
tonic vortex [18], i.e. the state containing a single vortex 
line, which have been recently observed in elongated ge¬ 
ometries pH Go]. It is worth to mention that also in 
toroidal systems, off-centered vortex lines (solitonic vor¬ 
tices) have been observed [8] in the density regions of the 
condensate. 

We have numerically solved the 3D GPE ([!]), by using 
the Newton method, in order to find steady states sup¬ 
porting stationary currents in toroidal condensates with 
typical experimental parameters. Fig. |3ja) displays data 
for a regime where, apart from winding number states, 
two types of solitary waves are possible: dark solitons 
and solitonic vortices. Our results show that 3D toroidal 
condensates support standing waves equivalent to those 
found in straight channels, and follow the same energy 
criterion for their bifurcation. At the scale of Fig. Up), 
the energies of dark solitons and solitonic vortices are in¬ 
distinguishable, and the main difference is their range of 
existence. While dark solitons there exist only for values 
of the angular momentum per particle around (q + l/2)h, 
solitonic vortices extend into the regions approaching 
winding numbers. This approach is mediated by the two 
different branches of Fig. |3|a) that are represented by 
means of lines with right and left pointing triangles, along 
with two characteristic density isocontours. The former 
branch is connecting with the winding number q = 0, 
and consists of vortices that produce flows rotating in 
the same direction as f2, while the latter connects with 
the winding number q — 1 and contains counterrotating 
vortices (antivortices). 



FIG. 3. Same as Fig. ^ for a 3D system, showing numeri¬ 
cal results of time-independent 3D GPE 0 for condensates 
containing 5.3 x 10 4 atoms of 23 Na in an isotropic harmonic 
trap with uj — x 100 Hz centered at R = 6a±. Two soli¬ 
tary waves can be excited in the absence of weak link (a), 
solitonic vortex (SV) and dark soliton (DS), whereas only a 
DS is possible at the center of a weak link (b) with depth 
U w i = 0.5 fi and width wq — 0.6 R. The insets represent 
density isocontours (at 5% of maximun density), coloured by 
phase, of SVs belonging to two separated branches containing 
either vortices or antivortices. The external diameter of the 
tori is 17 a±. 


It is noteworthy that hysteresis is also present in soli¬ 
tonic vortex states, since the two mentioned branches are 
separated by an energy barrier occupied by dark solitons 
(thick green line in Fig. |3ja)). This energy difference 
between solitonic states is responsible of their stability 
properties. As shown in Fig. [4j dark solitons decay into 
























5 



time (ms) 

FIG. 4. Evolution in real time, obtained by solving the 3D 
time-dependent GPE, of a dark soliton state, with the same 
parameters as in Fig. [3J after an initial Gaussian perturba¬ 
tion. Dimensionless azimuthal density an\(0) is plotted as 
a function of time. Density isocontours (at 5% of maximum 
density) are showed for times (from left to right) 0, 32 and 
160 ms. 

solitonic vortices, which are dynamically stable states. 
To show it, we have evolved in real time a dark soliton 
state made of 5.3 x 10 4 atoms of 23 Na in an isotropic har¬ 
monic trap with uj = 2tt x 100 Hz centered at R = 6a±. 
After adding Gaussian noise to the initial stationary soli¬ 
tonic state, we observe that it survives during, approxi¬ 
mately, the first 20 ms of the simulation, before decaying 
into another solitary wave, a solitonic vortex, and emit¬ 
ting phonons. The subsequent evolution shows the mo¬ 
tion of the emergent soliton with a different inertial mass 
than the initial state mm- 

As in the ID case, multidimensional solitons can exist 
only while their velocity around the torus is lower than 
the speed of sound. For isotropic harmonic trapping, the 
speed limit is given by the azimuthal speed of sound c(6) 
that can be calculated by [23] 

/Vo?) - m 

c( " ) = V 2m„«>) ■ (7) 

In Fig. |3ja), the thick dotted lines represent the values 
given by Eq. ([ 7 ]). The manifest agreement with numeri¬ 
cal solutions of GPE (lines with open triangles), permits 
Eq. 0 to be employed for predicting the range where 
solitons exist and hence for measuring the size of hystere¬ 
sis loops. 

A. Weak link 

The presence of the weak link introduces new effects in 
3D systems. The density depletion reduces not only the 
local speed of sound, as in ID systems, but also the local 
energy available for transverse excitations. As a result, 


some solitonic states that can be excited in regions with 
homogeneous density could not survive at the weak link. 
This is exactly the case represented in Fig. [3]^b), where 
solitonic vortices are excluded from the weak link region, 
and only dark solitons can be excited inside. As rotation 
starts, non-null currents are generated in winding-phase 
states. In analogy with the Meissner effect, the system 
produces a counter flow inside the weak link region in or¬ 
der to exclude the vorticity field induced by rotation. In 
a multidimensional system, a vortex-anti vortex pair pro¬ 
duces such a flow [24, [25] . The rotating weak link leads 
surface modes to enter into the torus by nucleating a vor¬ 
tex dipole along the radial direction nnm. This dipole 
yields a phase jump which is opposite to the background 
flow and has the required value for keeping a winding 
phase. The phase jump grows with increasing angular 
rotation, up to a critical point where the flow velocity in¬ 
side the weak link reaches the local speed of sound, and 
no further increase is allowed. At this point, which is 
also the end of the trajectory for solitons inside the weak 
link, the localized phase jump vanishes and a phase slip 
is produced. This effect is due to the dispersion curve 
of vortex dipoles [26 , which transform into phonons at 
low momentum p, that is to say, when E/p c. As a 
consequence, the size of hysteresis loops in 3D rings is de¬ 
termined by the range of existence of solitonic solutions, 
and this range is bounded by the local speed of sound 
c{6). 

Figure [5] shows the approach to the critical angular fre¬ 
quency followed by winding number states (upper graph 
in the figure) and solitons (lower) for a torus with the 
same parameters as before. At intermediate rotation 
rates, Cl = 0.5, a dark soliton has angular momentum 
per particle (L z )/Nh = 0.5, whereas the winding number 
q = 0 produces a small current, represented on the figure 
by an almost imperceptible gradient phase. However, as 
can be seen in the isocontours situated on the right at 
Fig. 0 at the rotation rate Cl = 1.87, close to the critical 
frequency, the differences between solitonic and winding 
numbers begin to disappear, and the converging process 
finishes with a phase slip to winding number q = 1, when 
the local speed of sound is reached. 


V. CONCLUSIONS 

We have studied the generation of persistent currents 
by means of rotating weak links in 3D Bose-Einstein con¬ 
densates confined in a torus. We demonstrate the exis¬ 
tence of persistent currents with non-quantized angular 
momentum supported by stationary states that include 
solitary waves. They are equivalent to the solitonic waves 
found in straight channels m, and follow the same en¬ 
ergy criterion for their bifurcation. In the presence of a 
rotating weak link, stationary solitons lying at the weak 
link give the lowest height of the energy barrier sepa¬ 
rating winding numbers and preventing phase slips. We 
show that the energy barrier vanishes when the angular 
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FIG. 5. Density isocontours (5% of maximum density) of 
winding number q — 0 (above) and solitonic (below) states, 
obtained by solving the time-independent 3D GPE [I] in 3D 
condensates with the same trap and rotating weak link as 
those of Fig. [3] Angular momentum per particle ( L z )/Nh 
is indicated in the inner part of the tori. The angular fre¬ 
quencies considered are Cl = 0.5 (left) and Cl = 1.87 (right), 
the latter near the critical frequency where states approach 
(. L z )/Nh — 0.18. Background graphs show the dimension¬ 
less density ani(0), after integrating along (z,r±). Colours 
represent the phase following the colour map of Fig. [3] 


frequency makes such solitons to move at the local speed 
of sound. At this rotation rate, the families of winding 
numbers and solitons converge, in such a way that time- 
independent states belonging to these families do not ex¬ 
ist for a faster rotation, and a phase slip is produced. 


Our results point out that the size of the hysteresis 
loops observed in Ref. [9] is characterized by the con¬ 
vergence process between winding numbers and solitons. 
This connection allows to identify, at the end of the hys¬ 
teresis loop, both stationary vortex lines (making soli¬ 
tonic states), and phonons (traveling at the speed of 
sound on winding number states). This fact indicates 
that Feynman and Landau criteria for the generation of 
phase slips meet at the connection point. 


However, the controversy about the different sizes of 
hysteresis loops given by Gross-Pitaevskii theory and ex¬ 
periments remains. As long as we are dealing with quan¬ 
tum mechanical systems, quantum tunneling processes 
can be considered in order to get phase slips [27]. It 
is not necessary for the system to wait for changing its 
angular momentum up to the configuration where the en¬ 
ergy barrier vanishes. These calculations are beyond the 
aim of the present work and will be reported elsewhere 

m- 
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